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Abstract. In 1969, Andreev and Lifshitz have conjectured the existence of a supersolid phase taking place 
at zero temperature between the quantum liquid and the solid. In this and a succeeding paper, we re-visit 
this issue for a few polarized electrons (spinless fermions) interacting via a U/r Coulomb repulsion on a 
two dimensional L x L square lattice with periodic boundary conditions and nearest neighbor hopping t. 
This paper is restricted to the magic number of particles N = 4 for which a square Wigner molecule is 
formed when U increases and to the size L — 6 suitable for exact numerical diagonalizations. When the 
Coulomb energy to kinetic energy ratio r s = UL/(2tV^N) reaches a value rf w 10, there is a level crossing 
between ground states of different momenta. Above r s , the mesoscopic crystallization proceeds through 
an intermediate regime (rf < r s < rY ~ 28) where unpaired fermions with a reduced Fermi energy 
co-exist with a strongly paired, nearly solid assembly. We suggest that this is the mesoscopic trace of the 
supersolid proposed by Andreev and Lifshitz. When a random substrate is included, the level crossing at 
rf is avoided and gives rise to a lower threshold rf(W) < rf where two usual approximations break down: 
the Wigner surmise for the distribution of the first energy excitation and the Hartree-Fock approximation 
for the ground state. 

PACS. 71.10.-w Theories and models of many-electron systems - 73. 21. La Quantum dots - 73.20.Qt 
Electron solids 



1 Introduction 

. A basic question in quantum many body theory is to know 
how one goes from independent particle motion towards 
collective motion when one decreases the density n s of a 
system of charged particles repelling each other via a U/r 
Coulomb repulsion. As introduced long ago by Wigner, 
the factor r s , defined as the radius r of the volume enclos- 
ing a single particle in units of the Bohr radius a b , governs 
this crossover. Quantum mechanical effects are important 
when r s is small, and become more and more negligible 
when r s becomes large. Since r s oc l/y/n s in two dimen- 
sions, the quantum limit of a Fermi liquid is obtained at 
large densities n s . In the dilute limit, the quantum effects 
disappear and the charges crystallize, forming a Wigner 
solid of minimum electrostatic energy. A calculation of the 
electrostatic energy of different crystalline arrays shows 
that the hexagonal array minimizes the energy in the 2d 
continuum. However, in a square lattice model with peri- 
odic boundary conditions (BCs) and for a sufficient filling 
factor, the symmetry of the Wigner solid is restricted and 
the formation of a square crystalline array is favored. 

Before studying this quantum-classical crossover in a 
mesoscopic lattice model, let us mention that it is usually 
assumed |1I2| that a single liquid-solid transition takes 
place at r s ~ 37 in the continuous 2d thermodynamic 



limit. This single transition was obtained by fixed node 
quantum Monte Carlo calculations |1I3| , allowing to study 
a few hundreds of electrons and to vary their number for 
estimating the finite size effects. When the spin degrees 
of freedom are included, the existence of an intermedi- 
ate polarized liquid phase separating the unpolarized liq- 
uid and the Wigner solid is still debated 0]. Though the 
quantum Monte Carlo calculations have the advantage to 
allow the study of a relatively large number of particles, 
they have the well-known "sign problem" . This leads to 
fixed node approximations which are made to avoid the 
negative weights that would be generated otherwise by 
antisymmetric states, and gives only an upper bound to 
the exact ground state energy. In Refs. jll'ij for instance, 
two nodal structures have been considered, given by two 
Slater- Jastrow wave functions adapted to describe the 
weak coupling Fermi liquid (nodal structure of a Slater de- 
terminant of plane waves) and the strong coupling Wigner 
solid (nodal structure of a Slater determinant of localized 
site orbitals). In these works, the existence of a single tran- 
sition separating the Fermi liquid from the Wigner solid 
is a consequence of the assumed nodal structures. Very 
recently, the ground state energy of the two-dimensional 
uniform electron gas has been calculated [H] with a fixed 
node diffusion Monte Carlo method, including backflow 
correlations. The backflow method allows to partly re- 
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lax the previous nodal constraints. The backflow nodes 
give smaller energies that the previous nodes, supporting 
the existence of an intermediate polarized liquid phase for 
26 < r s < 35. 

1.1 2d metal and related issues 

The motivation to re- visit nowadays charge crystallization 
in 2d Coulomb systems is fourfold. Firstly, it becomes pos- 
sible to create 2d gases of charges in high quality field ef- 
fect devices and to decrease by a gate the carrier density n s 
for obtaining a large factor r s . Doped semi-conductors 
(Si-Mosfet, Ga-As heterostructures, Si-Ge quantum wells) 
can be now used to study how one goes in two dimensions 
from a Fermi liquid towards a Wigner crystal. Secondly, 
the direct observation of the Wigner crystal being diffi- 
cult, one can nevertheless measure the conductance of the 
dilute 2d electron gas at different densities as a function of 
the temperature, of the bias voltage, of a parallel magnetic 
field, etc. Remarkably, those transport measurements, first 
done [H] by Kravchenko et al in high quality Si-MOSFETs, 
show the existence of an unexpected metal-insulator tran- 
sition (MIT) when the gate voltage varies and the dilute 
limit is reached. A similar MIT was later observed using 
a hole gas in Ga-As heterostructures and in Si-Ge quan- 
tum wells. This observation of a low temperature metallic 
behavior for typically 3 ... 6 < r s < 9 ... 30 (the largest 
ratios characterizing the cleaner samples) raises the ques- 
tion of a possible intermediate phase, which should be 
neither a Fermi glass of localized particles (Anderson in- 
sulator), nor a correlated Wigner solid (pinned insulating 
crystal). Quite recently, magnetotransport measurements 
in low density 2d hole gas in GaAs quantum wells have 
been interpreted |7I8) in terms of an unknown metallic 
phase coexisting with a Fermi liquid phase, for estimated 
values 12 < r s < 18. Local compressibility measurements 
by Hani et al also point to a two-phase coexistence pic- 
ture for intermediate values of r s in GaAs HJ. Thirdly, 
the formation of a mesoscopic Wigner molecule can be 
also nowadays studied using a quantum dot^Oj with a 
few electrons or a few ions [IJ trapped by electric and 
magnetic fields. Increasing the size of the trap yields 
a crossover from independent-particle towards collective 
motion. Quantum dots with a few electrons are among 
possible candidates for providing the Qubits of a future 
quantum computer |13) . Lastly, an unexplained intermedi- 
ate regime was numerically observed |14I15I16*] in studying 
the persistent currents carried by the ground state and the 
low energy excitations of mesoscopic disordered clusters. 
Both experiments and numerics give unexplained low en- 
ergy behaviors for similar intermediate values of the ratio 



1.2 Andreev-Lifshitz supersolid 

In this study, we consider fully polarized electrons and 
ignore possible magnetic transitions. For avoiding un- 
controlled assumptions, we consider a system which is 



small enough to allow exact numerical diagonalization, 
but where the lattice effects and the finite size correc- 
tions are important. Those two effects have been stud- 
ied in details in Ref. ^Zj for two polarized electrons. Our 
main goal is to characterize as precisely as possible the 
ground state for intermediate ratios r s , in order to see if 
a small lattice model does not exhibit the mesoscopic sig- 
nature of an intermediate phase separating the solid from 
the liquid, where the solid and the fluid would coexist. 
Such a vacancy-solid phase was indeed suggested ^Hj by 
Andreev and Lifshitz if the zero point motions of certain 
defects become sufficient to form waves propagating inside 
the solid. This Andreev-Lifshitz supersolid was first pro- 
posed for three dimensional quantum solids made of He 
atoms. Castaing and Nozieres have later considered [B?] 
such a possibility for spin polarized He 3 . The statistics of 
the defects depend on their nature. For simple vacancies 
in the crystal, their statistics is given by the statistics of 
the particles out of which the solid is made. If the defects 
are bosons, they may form a condensate, giving rise to a 
superfluid coexisting with the solid. This supersolid phase 
is discussed in certain bosonic models [201 • If the defects 
are fermions, they may form a Fermi liquid |21| coexist- 
ing with the solid, such that the system is neither a solid, 
nor a liquid. Two kinds of motion are possible in it; one 
possesses the properties of motion in an elastic solid, the 
second possesses the properties of motion in a liquid. 

If one considers the quantum melting of the solid from 
the dilute limit (large r s ), the nature of the relevant de- 
fects is not an easy question. One can imagine a particle 
being put into an intersticial site of the Wigner lattice, cre- 
ating a vacancy-intersticial pair at a certain electrostatic 
cost SU. A model assuming such defects has been recently 
proposed 22 for describing addition spectra in quantum 
dots. Classically, this vacancy-intersticial pair remains lo- 
calized. But quantum tunneling may lead to derealization 
of the defects and to the appearance of a band of defects of 
finite width Bd which increases when r s decreases. When 
Bd exceeds SU, one can imagine two possibilities: either 
the total quantum melting |5] of the Wigner crystal (sim- 
ple solid- liquid transition), or a partial melting leading to 
the persistence of a floppy crystal with delocalized defects. 
If a delocalized defect appears in the quantum crystal, 
the crystal remains perfectly periodic, but the number of 
crystal lattice sites becomes smaller than the total num- 
ber of particles. This is the supersolid scenario proposed 
by Andreev-Lifshitz for He physics. 

If one considers charge crystallization from the other 
limit, where the density n s is large (small r s ), one can 
argue that the interaction will create correlated pairs of 
particles near the Fermi surface, but will not reorganize 
the one particle states well below the Fermi surface. Such 
a possibility has been proposed by Bouchaud et al (231 
for liquid He 3 . Moreover, they have developed a varia- 
tional approach, based on a fixed number of fermions BCS 
wave function, having a different nodal structure than the 
Jastrow-Slater nodal structures considered in Ref. pQ. In 
this picture, the system is thought as made of unpaired 
fermions with a reduced Fermi energy, co-existing with a 
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strongly paired, nearly solid assembly. Furthermore, it was 
stressed that the supersolid - first introduced by Andreev 
and Lifshitz - would be a good candidate to describe this 
new phase. To the concept of a crystal with a reduced 
number of crystal lattice sites, as discussed by Andreev 
and Lifshitz from the solid limit, corresponds the concept 
of unpaired fermions with a reduced Fermi energy, as dis- 
cussed by Bouchaud et al from the liquid limit. 

Our purpose is to study if a supersolid regime is rele- 
vant for describing a low density 2d electron (or hole) gas. 
This numerical work being restricted to very small system 
sizes, the study is limited to merely study if one can de- 
tect the mesoscopic trace of a possible supersolid. One of 
our motivations comes from the observation that the 
low energy levels do not obey Wigner-Dyson statistics for 
disordered clusters at intermediate ratios r s . This will be 
again emphasized at the end of this work for the statistics 
of the first energy excitation. This suggests the existence 
of low energy collective excitations in the clean limit. For 
this reason, we have studied clusters without random sub- 
strate and we have observed at intermediate r s a floppy 
correlated solid coexisting with a liquid of unpaired par- 
ticles. This conclusion is supported by a study of the pro- 
jection of the ground state (GS) onto a combination of 
Slater determinants (SDs) built out from plane waves and 
from site orbitals. The plane wave SDs are given by the 
low energy levels of same total momentum K as the inter- 
mediate GS, and correspond to unpaired fermions with a 
reduced Fermi energy. The site SDs describe the Wigner 
solid molecule and its small fluctuations. Since the GS is 
given by the combination of unpaired fermions and of a 
floppy Wigner molecule for rf « 9.3 < r s < rY ~ 28 
in the studied system, we suggest that this is the meso- 
scopic trace of the supersolid discussed in Refs. and 
[231 - The study of the GS response to various perturba- 
tions (Aharonov-Bohm flux, pinning well) and of the dis- 
tributions of the different inter-particle spacings allows us 
to give a few remarkable properties of the intermediate 
regime. Eventually, we consider the effect of disorder and 
give further evidence of the existence of a first thresh- 
old rf(W) where strong correlation effects occur without 
yielding a full crystallization: (i) the breakdown of the 
Hartree-Fock approximation for the ground state and (ii) 
the breakdown of Wigner-Dyson level repulsion for the 
first excitation. 



erators of a spinless fermion at the site i; n\ = c\c{ is 
the occupation number at the site labeled by the vector 
i = (i x ,i y ). The vector ry is defined as the shortest vec- 
tor going from the site i to the site j in a square lattice 
with periodic BCs (r x and r y < L/2). t — h 2 /(2ma 2 ) is 
the hopping term, a the lattice spacing, v\ the site po- 
tentials which are randomly distributed in the interval 
[-W/2,W/2] and U = e 2 /(ea) the Coulomb interaction 
between two fermions separated by a in a medium of di- 
electric constant e. 

This work is restricted to a detailed study of the case 
N = 4 and L — 6, corresponding to a filling factor v = 
N/L 2 = 1/9. N — 4 is a 'magic' number for which at 
large values of U, the 3x3 square Wigner molecule is 
commensurate with the imposed 6x6 square lattice. The 
r. factor, defined in the continuum as 



(2) 



for a carrier density n s and a Bohr radius as — h 2 e/(me 2 ), 
becomes in a lattice model 



U 



2t 



(3) 



since h 2 /(2ma 2 ) — > t, e 2 /(ea) — ► U. and n s — v/a 2 . 

A L x L continuous 2d torus having infinitely more 
degrees of freedom than a mere 6x6 lattice, one can- 
not compare the obtained lattice behaviors to those ob- 
tained assuming a continuous space, as in Ref. pQ, without 
further investigations. Nevertheless, appropriately defined 
observables should only depend on the value of the di- 
mensionless ratio r s , up to certain finite size corrections. 
This has been checked ^7] for two polarized electrons on 
a square lattice, when the fluctuation Ar of the distance 
r between the two particles is larger than the lattice spac- 
ing a. Calculating Ar in powers of t/U, one finds j^TJ that 
a correlated lattice regime takes place when r s exceeds a 
threshold value r* w 100 when N — 2 and L = 6. Below 
r*,r s oc UL/t is the relevant scaling variable, up to certain 
finite size corrections of order 1/L 2 . Above r*, one has a 
lattice regime where r s is not a relevant scaling variable. If 
the threshold value r* does not vary very much when one 
goes from N = 2 to N = 4, one has a chance to observe 
a four particle Wigner molecule free of important lattice 
effects as far as r s < r* w 100 when L = 6. 



2 Lattice model 

We consider fully polarized electrons (i.e. spinless 
fermions), having symmetric spin wave functions and an- 
tisymmetric orbital wave functions, free to move in an 
Lx L lattice with periodic BCs, and interacting via aU/r 
Coulomb repulsion. The Hamiltonian reads 

<ij> i t*j 1 iji 

where i, j label the lattice sites, means i nearest 

neighbor to j, c\, c\ are the creation, annihilation op- 



3 The non disordered lattice 

When there is no disorder (W = 0), k = (k x ,k y ) being 
the one particle momentum, it is more convenient to write 
H using the Fourier transforms of the creation and anni- 
hilation operators. One has the relations 

cj = ^£^ lk - j , (4) 

k 

and 

d k = jJ2 c ^- ( 5 ) 
j 
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Fig. 1. Low kinetic energy plane wave SDs in momentum space 
of coordinates (k x ,k y ): 1 GS \Ko{/3)) with K/0 (upper left), 
together with 3 plane wave SDs of K = and of low energies: 
one of the 4 \Ki((3)} of energy — 12t (upper right), and two 
of the 16 ^4(7)) of energy — 9t (lower right and left) directly 
coupled by the two body interaction to the \Ki(/3)) shown in 
the upper right figure. A circle means that the state (k x , k y ) is 
occupied, its energy being indicated in units of t. 
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Fig. 2. Low Coulomb energy site SDs in real space of co- 
ordinates (x,y): one \S) with its 9 possible centers of mass 
(upper left) , two IP2} (upper right), two IP3} (lower left) and 
two \DS) (lower right), which give by successive translations 
one |S(K = 0)}, two |P 2 (K = 0, J)), two |Pi(K — 0, J)} and 
two out of four \DS(J£ = 0, J)) respectively. The two others 
\DS(K = 0, J)} are obtained by successive translations of the 
square deformed at the opposite corner. 



which yield 

ff=^4d k£ (k)+ Y, ^(q^L+q^L-q^rfk, (6) 
k q,ki,k 2 

where 

e(k) = — 2* (cos k x + cos k y ) (7) 

and 

In the eigenbasis of the non interacting system (eigen- 
vectors dj^dj^dk |0), |0) being the vacuum state), the 
Hamiltonian matrix is block diagonal, each block being 
characterized by the same conserved total momentum 
K = k^. Only the non interacting states having in 

common two ks out of four can be coupled by the interac- 
tion inside a K sub-block. Therefore, each K sub-block is 
a sparse matrix which can be exactly diagonalized using 
the Lanczos algorithm. 

3.1 The free Fermi limit 

When U = 0, the states are then Njj plane wave SDs 
d^ljlA. I°>- n b = M\/{N\{M - N)\) = 58905 for 
M — L 2 — 36 and N = 4. The low energy levels without 
interaction are by increasing energies: 

- 4 GSs \K Q {P)) of energy E (U = 0) = -13* and of 
momenta Kq 7^ 0. 



— 25 first excitations of energy Ei(U = 0) = —12*, 

— 64 second excitations of energy Ei(U = 0) = — 11* and 
of momenta K2 7^ 0. 

— 180 third excitations of energy E^(U = 0) = —10* and 
of momenta K3 7^ 0. 

— 384 fourth excitations of energy E^U = 0) = —9*. 

If one considers the low energy states of total momen- 
tum K = 0, some of them being shown in Fig.^ one finds 
by increasing energy: 

— 4 SDs \Kx{fi)) ((3 = 1,...,4) of energy -12*, corre- 
sponding to a particle at an energy —At with ki = 
(0,0), two particles at an energy —3* and a fourth 
particle of energy —2*; plus a single SD \Kx(0)) with 
4 particles of energy —3*. 

-16 SDs ^4(7)) (7 = 1,...,16) of energy -9* 
given by 8 SDs where the particles have energies 
—4*, —3*, —2*, 0* respectively and by 8 other SDs where 
the particles have energies —3*, —3*, —2*, — * respec- 
tively. Note that the ^4(7)) are directly coupled to 
the \Ki(f3)) by the pairwise interaction. 

3.2 The correlated lattice limit 

When * = 0, the translational invariance is broken and 
the states are Njj Slater determinants Cjcjc^cj |0) built 
out from the site orbitals. The configurations ijkl corre- 
spond to the Njj different patterns characterizing 4 dif- 
ferent sites of the Lx L square lattice. The configurations 
of low electrostatic energy are respectively: 



G. Katomeris et al: Andreev-Lifshitz supersolid I 



5 



— 9 square configurations \S) of side b = 3 and of energy 
E (t = 0) re 1.8017, 

— 36 parallelograms \P\) of sides (3, \/l0) and of energy 
~ 1-85C/, _ _ 

— 36 other parallelograms \P^) of sides (vlO, vlO) and 
of energy ps 1.97/7, 

— 144 deformed squares \DS) obtained by moving a sin- 
gle site of a square \S) by one lattice spacing and of 
energy pa 2U . 

Some of those low energy site SDs are shown in Fig. El 
When an infinitesimal hopping term t is included, one 
must delocalize the site SDs in order to restore transla- 
tional invariance and to have eigenstates of given quan- 
tized total momenta K. For instance, the 9 squares \S) 
give 9 eigenstates of momentum K 



1 L 



(9) 




where 



Tj \S) = c\ 



|o>. (io) 



The possible momenta for the IS'(K)) are given by 
(K x ,K y ) = 27r(j) x , Py )/(L/2). 

If we consider the low energy states of total momentum 
K = when r s — > oc, one finds by increasing energy the 
following delocalized site SDs: 

— 1 delocalized square \S(K = 0)), 

— 2 delocalized parallelograms |Pi(K = 0, J)) (J = 1, 2) 
obtained from the 36 \Pi), 

— 2 delocalized parallelograms = 0, J)) (J = 1, 2) 
obtained from the 36 IP2), 

— 4 delocalized deformed squares \DS(K = 0, J)) ob- 
tained from the 144 \DS). 



3.3 Level crossing at r l s 



and charge crystallization at 



The low energy part of the spectrum is shown in Fig. ^ 
as a function of r s . If we follow the 4 GSs of energy — 13i 
at r s = (Ko 7^ 0), one can see a first level crossing at 
rf w 10 with a non degenerate level (Ko = 0) which be- 
comes the GS above rf, followed by two other crossings 
with two other sets of 4 levels with Kr 0. When r s 
is large, 9 levels coming from Ei(r s = 0) have a smaller 
energy than the 4 levels coming from Eo(r s = 0). Since 
the degeneracies are (9, 36, 36) when t = 0, these 9 states 
give the 9 square molecules \S > when r s — > 00. The de- 
generacies ordered by increasing energy become (1, 4, 4, 4) 
instead of (4, 25, 64) for r s = 0. 

To describe large r s , one can use degenerate pertur- 
bation theory and study how the degeneracy of the 9 \S) 
is removed by terms oc t/U oc r.T 1 - The centers of mass 
R of the 9 \S) are located on the periodic 3x3 square 
lattice sketched in Fig. |3 For large r s , one has a single 
rigid molecule free to move on this restricted lattice, with 



Fig. 3. As a function of r s , low energy part of the spectrum 
exhibiting a level crossing at rf . Inset: two first level spacings 
Ai/t (dashed) and Aijt (dotted) and the perturbative result 
3Dr~ 3 (dot-dashed). 



a hopping term T oc trj 3 and a quantized momentum 
K = 2tt/(L/2)(p x ,p v ). Taking into account also the cor- 
rections to the diagonal matrix elements, one obtains for 
the 9 first energies Eo(K) in the limit r s — > 00 

= E D - 2j(cos K X (I)+ cos K y (I)). (11) 

E D = Ar s + B/r s + C/r 3 s and T/t = D/r 3 s (A ps 2.13, 
B ps -70.81, C ps -18763 and D ps 3464). E D comes from 
the small vibrations of the rigid molecule while 8T is the 
band width of its zero point fluctuations. The degeneracies 
are 1, 4, 4 respectively. 

Four observations can be drawn from this t/U expan- 
sion. 

— The ground state must exhibit a level crossing since the 
total momentum K = when r s — > 00 (lowest quan- 
tized kinetic energy for the center of mass of a rigid 
square molecule) while K ^ when r s — > (incom- 
plete filling of the Fermi shell — 3i). Is this GS level 
crossing a general feature? As explained in Ref. 
K / at large r s for N = 3 and L = 6, and there 
is no GS level crossing while there is one if = 3 
and L = 8. When the spins are included, the GS level 
crossing disappears for N — 4 and L = 6. K = at 
large r s for any L x L square lattice with a filling fac- 
tor v = 1/9. For L = N = 9, the Fermi shell -2t is 
totally filled, K = at r s — and momentum conser- 
vation does not yield a GS level crossing, in contrast 
to the case L = 12 and N = 16 where the Fermi shell 
is incompletely filled. As we sec, the existence of a GS 
level crossing depends on L and N and may not have 
a particular significance. In this work, we have studied 
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the true GS, taking the subspaces of K ^ below rf , 
the subspace of K = above rf . One could have pre- 
ferred to study the GS inside the subspace of K = 
for all the values of r s , to find that the onset of cor- 
relation effects which we observe at rf as we will see 
later should occur at a possibly smaller threshold in 
the K = subspace. 

— In the inset of Fig- El one can see that the t/U expan- 
sion gives an accurate description of the 9 first energies 
above a relatively large value r s sa 100. As explained 
in Ref. J7|, this t/U expansion is characteristic of a 
correlated lattice regime where the fluctuations of the 
charges around the equilibrium Wigner lattice sites are 
strongly restricted by the imposed lattice. This lat- 
tice expansion has to be distinguished from the large 
r s expansion of a continuous model, where the oscilla- 
tory motion of the electrons around the Wigner crystal 
equilibrium positions gives [25] for the GS energy an 
expansion in powers of 1/^/rJ. 

— Though the t/U lattice expansion ceases to be accu- 
rate below r s « 100, the 9 low energy states begin to 
have the structure of the spectrum of a single mas- 
sive particle in a 3 x 3 lattice (two equal energy spac- 
ings Ai — A2 characterizing the 3 first sets of states 
with degeneracies 1,4,4 respectively) at a lower value 
ry ~ 28. This structure means that the system re- 
mains essentially a rigid square molecule with its 9 
quantized modes for the motion of the center of mass 
down to rY ■ To create a defect in this square molecule 
costs a high energy available in the 10 t/l excitation 
only. This is why we identify rf as the threshold value 
for the mesoscopic crystallization, above a first thresh- 
old rf and below a higher threshold r s w 100 where 
the lattice t/U expansion becomes valid. 



3.4 Truncated site basis 

The N H site SDs c|c]44 |0) correspond to the N H dif- 
ferent patterns characterizing 4 different sites ijkl of the 
L x L square lattice. If we order those configurations 
by the smallest distance d between two sites, Nd denot- 
ing the number of configurations with inter-site spacings 
larger than d, one has ATj = 27225, A^ = 9837, N 2 = 
2709, = 81 configurations having a smallest inter-site 
spacing > d, out of Nh = 58905 configurations. The two 



thresholds rf and 



can be also detected if one calcu- 



lates the GS energy Eo(Nd) of the truncated Hamiltonian 
written using the site SDs basis restricted to Nd site SDs 
and if we consider the error AE (Nd) = (E (Nd) — E )/t 
made using this truncation for having the exact GS energy 
Eq. As shown in Fig. [I] the error AEq{Ni) becomes small 
above rf , while the error AE (N^) = Eo(t = 0) — E for 
all values of r s . The error AEo(N2) has a very interest- 
ing behavior. As r s increases, AE (N 2 ) first decreases up 
to r s ~ rf , then exhibits a very remarkable plateau for 
rf < r s < r^ , taking a value w it independently of r s , 
before decreasing as AEq(N/^) above rf ■ This plateau 
suggests that the GS for intermediate r s is composed of a 
floppy molecule which can be projected onto the N 2 site 




Fig. 4. 

r s : d = 
(□). 



Errors AE (N d ) 
1 (thick line), d = 



= (Eo(Nd) — Eo)/t as a function of 
V2 (dotted line), d = 2 (o), d = 




Fig. 5. GS Projections P ° (o), (□), P 2 (♦), P 3 (x) and 
Po (*) onto plane wave SDs of low energy when r s — > 



SDs adapted to describe it, plus an unpaired fermion of 
kinetic energy « — 3i which is not included in this trun- 
cated subspace since it is delocalized. Very remarkably, 
this energy turns out to be the energy of a particle at the 
Fermi surface of the non interacting system. 



3.5 GS projections onto plane waves 

To understand further the nature of the GS between rf 
and rf , we have projected the GS wave functions |$b( r s)) 
over the low energy plane wave SDs appropriate to de- 
scribe unpaired fermions. As shown in Fig.0 below rf , a 
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Fig. 6. GS projection P^(r s ) (■), P^n) (*), P^{r a ) (0) 
and P^o(r s ) (<l) onto the first K = delocalized site SDs of 
low energy when r s — > oo. 

K^O GS has a large projection Pq (r s ) over the ?7 = GS 
of same K and begins to have a smaller projection P 2 (r s ) 
over the second excitations of the non interacting system 
of same K. Above rf, the non degenerate GS with K = 
has of course no projection onto the plane wave SDs of 
K/0, but has a large projection 

n 1 (^) = El(^(r s )|ifi(/3))| 2 (12) 

which is equally distributed over the 4 low energy states 
\Ki(f3)) of K = and a non negligible projection 

16 

^o(^) = El^o(r s )|X4(7)>| 2 (13) 

7=1 

over the 16 previously defined ^4(7)) of K = which are 
directly coupled by the interaction to the \K\ (/?))• 

One concludes that a significant part of the system 
remains an excited liquid above rf , described by a large 
projection Pq = Pq + Pq over a few combinations of low 
energy unpaired fermions. Due to the GS level crossing, 
the intermediate GS has to be described from the K = 
Fermi sea and not from the K ^ Fermi sea. Since this 
projection is only partial, only a part of the system is 
made of unpaired fermions, in agreement with the concept 
proposed by Bouchaud et al of a reduced Fermi energy, 
which decreases as r s increases. 

3.6 GS projections onto site orbitals 

We study now the GS projection over the low energy 
site orbitals shown in Fig. [3 which become the eigen- 
states when r s — > 00. More precisely, we consider the 



first delocalized site SDs having a delocalized center of 
mass and the same momentum K than the GS: The de- 
localized square IS'(K)), the 2 delocalized parallelograms 
|Pi(K,J)), the 2 delocalized parallelograms \P2CK, J)), 
and the 4 delocalized deformed squares \DS(K 7 J)). Fig. 
shows the behaviors of the GS projections 

P^(r s ) = \{^ (r s )\S(K))\ 2 (14) 
2 

^(^)-Ei(^^)i p i( K ' J ))i 2 ( 15 ) 
j=i 

2 

^(^) = El^^)l P 2(K,J))| 2 (16) 

j=i 

^)H(^(r s )|PS(K,J))| 2 , (17) 

where K is the GS momentum (K 7^ below rf and 
K = above rf). 

While aK/0 GS has negligible projections over the 
low energy site SDs of same K below rf , there is an im- 
portant contribution above rf of the deformed squares, 
of the square and of the parallelograms 1 of K = 0. As 
r s increases, the GS projection P^(r s ) over the square 
molecule of K = goes to 1. The GS projection P^ D (r s ) 
over the 4 deformed squares of K = is the main pro- 
jection below r s m r^ , a threshold value where the GS 
projection P^(r s ) over the parallelograms 1 is maximum. 
One concludes that above rf , the missing part of the 
system, which is not described by the low energy un- 
paired fermions of the previous section, is a floppy Wigner 
molecule, mainly made of deformed squares below r^ and 
of squares above r^ 7 . 

3.7 GS projections onto a combined basis of plane 
waves and site orbitals 

The site SDs and plane wave SDs are not orthonormal. Af- 
ter re-orthonormalization, the total GS projection Pt(r s ) 
over the subspace spanned by the 20 plane wave SDs and 
the 9 delocalized site SDs of low energy and of momentum 
K = is given in Fig. Q together with the GS projection 
Po( r s) = Po( r s) + Po( r s) over the 20 plane wave SDs and 
Px>( r s) = S/=o Poo( r s) over the 9 delocalized site SDs of 
momentum K = 0. One can see than more than 95/100 
of the intermediate GS is located inside this combined 
subspace, suitable to describe a floppy solid co-existing 
with low energy unpaired fermions. This demonstrates the 
Andreev-Lifshitz conjecture for the considered mesoscopic 
lattice model. 

Let us point out that those exact results raise serious 
objections about the validity of "exact" studies of a few 
particles in the continuum, where the infinite Hilbcrt space 
is truncated to the finite basis made of the low energy 
plane wave SDs only. For intermediate r s , our exact results 
show the necessity to combine plane wave SDs and site 
SDs before truncating. 
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Fig. 7. GS projection Po(r a ) (*) over the 20 K = plane wave 
SDs of low energy, Poo{r s ) (□) over the 9 K = delocalized 
site SDs and Pt(r s ) (•) over combined plane wave and site SD 
re-orthonormalized basis. 
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Fig. 9. Disorder average of the relative fluctuations (ud(i)} 
of the 6 inter-particle spacings as a function of r s , using the 
same symbols as in Fig. |H| Characteristic behaviors u r — 
0.28— 0.01r s (dotted-dashed), u T oc 6/r s (solid) and u r oc 26/r^ 
(dashed) . Note the three spacings having an almost r s indepen- 
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Fig. 8. Disorder average of the 6 mean inter-particle spacings 
(d(p)) as a function of r B . 

3.8 Inter-particle spacings 

To understand the nature of the intermediate GS, we 
study the distribution of the different inter-particle spac- 
ings. For the site SDs Cjcjc^cJ |0), one defines the 6 spac- 
ings e?ijki(l) < rfyki(2) < . . . < rfyki(6) of each configu- 
ration ijkl ordered by increasing values. The n th moment 
d n (p) of the p th GS inter-particle spacing at r s is given 
by: 

N H 

(d n ( P )) = J2 ^kiWI (*o(r.)| 4c\clc\ |0) I 2 . (18) 
ijki=i 



Very weak random potentials are included to get rid of 
the symmetries of the 6x6 lattice. After average over 
an ensemble of random configurations, we show in Fig. 
El how the 6 mean GS inter-particle spacings (d p ) vary 
as a function of r s for a value W = 0.1 of the disorder 
strength. When r s — > 00, the 3x3 Wigner molecule gives 
d\ = di = e?3 = d± = 3 and d§ = d§ = 3\/2- As r s decays, 
one can see that one of the largest spacings out of two 
and two of the smallest spacings out of four remain close 
their asymptotic values, in contrast to the others. This 
shows us that one has for intermediate r s a floppy solid 
made of three particles, while the fourth particle remains 
delocalized. A similar conclusion was drawn from a study 
of the case N — 3 and L — 6 in Ref. [3^, where the 
intermediate GS was shown to be a floppy two particle 
molecule co-existing with a third delocalized particle. 
The behaviors of the relative fluctuations 

-w-i/ilF (19) 

of the 6 inter-particle spacing d(p) are given in Fig. [5] 
When r s — > 00, the fluctuations of the square molecule 
can be calculated using the t/U lattice expansion. At first 
order, one can move only a single particle, which modifies 
three inter-particle spacings out of six. The fluctuations of 
the three remaining spacings is obtained by moving two 
particles, which requires to go to the second order. This 
explains the three r~ x decays and the three r~ 2 decays 
characterizing the correlated lattice regime. The behaviors 
in the intermediate regime are remarkable: 

— The relative fluctuations of three spacings out of six 
decay as r s increases, as one can expect if a floppy 3 
particle molecule becomes more rigid as r s increases. 
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Fig. 10. Functions a (r s )/t (dotted) and 6(r s )/t (filled squares) 
characterizing the GS response to an infinitesimal flux; inset : 
b(r s )/t (dashed) and 8D/9r s " 3 (dotted). 



The three others spacings have relative fluctuations 
which are nearly independent of r s between rf and 
, as one can expect if the 4 th particle remains delo- 
calized. 

Notably, the smallest spacing has a relative fluctuation 
which varies as 0.28— 0.01r s for r s < rf (weak coupling 
regime), which is almost independent of r s for rf < 



< r^ (intermediate regime) before decreasing as 
6/r s above r^ (correlated lattice regime). 

The behaviors of the different inter-particle spacings 
are consistent with the Andreev-Lifshitz conjecture for the 
considered mesoscopic lattice model. Notably, the number 
of crystal lattice sites is indeed smaller than the total num- 
ber of particles. 



3.9 GS response to small perturbations 

We study now the GS response to small perturbations, 
when the site potentials are non random (W = 0), which 
gives us other signatures of the intermediate regime. 




Fig. 11. GS density D p (r s ) at the pinning site p with V p /t 
-0.01 (dotted), -0.05 (dashed) and -0.1 (dot-dashed). 



crossover between two regimes around (inset of Fig. 

EE). 



3.9.2 Single pinning well 

The second perturbation consists in introducing a weak 
negative potential V p at a single lattice site p. The GS 
density 

D p (r s ) = (Mr s )\clc p \%(r s )) (20) 

at the site p is shown in Fig. ^2 If Vp = 0, D p (r s = 0) = 
1/9. A weak negative value of V p yields a larger value for 
D p (r s = 0). When one turns on the interaction, D p first 
increases and drops at rf , where the interacting GS begins 
to have a weaker response to a weak pinning well than the 
non interacting GS. When r s is large, D p increases again 
and the Wigner molecule is pinned. This surprisingly weak 
response for intermediate r s suggests that the system may 
very weakly respond to the presence of weak impurities. 
Let us underline that this is precisely for those values of 
r s that the new 2d metal has been observed in 2d field 
effect devices UJ. 



3.9.1 Aharonov-Bohm flux 

The first one consists in piercing the 2d torus by an in- 
finitesimal positive flux 4> (periodic BCs along the y di- 
rection, t — > texp(i(f>/L) for hopping along the x-direction 
only, <j> = 7T corresponding to anti-periodic BCs). The coef- 
ficients a(r s ) and b(r s ) (Kohn curvature) of the expansion 
E {r s ,4>) « E {r s ,Q) + a(r s )<f> + b(r s )(j) 2 /2 are given in 
Fig. ED When r s = 0, <j> removes the fourfold degener- 
acy of Eq, a = —V3t/6 and b — 7t/36. When r s is large, 
the substitution K X {I) -> K X (I) + 20/3 in Eg. ITT1 gives 
a = and b s» 8Dtr~ 3 /9. An infinitesimal positive flux 
4> gives rise to a persistent current I x = —dE /d<f) = —a 
when r s < rf while the GS curvature b exhibits a smooth 



4 Lattice model with random potentials 

We extend our study of the disorder by increasing W from 
the previously considered weak value (W = 0.1) up to 
larger values (W — > 5) which are too small for having 
exponential localization of the one particle states on a 
scale L = 6, but sufficient for having one particle diffusion. 

4.1 Inter-particle spacings 

In Fig. El one can compare the mean smallest and largest 
inter-particle spacings. For W — 0.1, the effect of the 
level crossing at rf is still visible when one follows the 
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Fig. 12. Disorder average of the mean smallest and largest 
inter-particle spacings for W = 0.1 (filled symbols) and 5 
(empty symbols) as a function of r s . 
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Fig. 14. Ensemble average of the first energy spacing (AE) = 
(Ex - £t,} as a function of r s for W = 0.2 (A), 0.5 (o), 1 (o), 
3 (V) and 5 (□). Inset: crossover values rf(W) as a function 
of W. 
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Fig. 13. Disorder average of the relative fluctuations of the 
smallest and largest inter-particle spacings as a function of r a 
for W = 0.1 and 5 (same symbols as in Fig. 1121 . 



largest spacing, but becomes very smooth for the small- 
est spacing. The jump associated to rf is totally smeared 
for W = 5. One can also see that the smallest spacing is 
closer to its asymptotic value 3 as one increases W, show- 
ing that disorder favors the formation of a correlated glass. 
In contrast, the disorder defavors the formation of a per- 
fect crystalline array, the largest spacing requiring a larger 
U to reach its asymptotic value 3^/2. Similar conclusion 
can be drawn from Fig. 1131 where the corresponding rela- 
tive fluctuations of the smallest and largest spacings are 
given. 



4.2 First energy excitation 

Fig. El gives the average over the disorder of the first en- 
ergy spacing (AE) = (Ex — E ) as a function of r s for 
different values of W. One can see a weak coupling regime 
where (AE) does not depend on r s , followed by an increase 
of (AE) as a function of r s . The threshold coupling be- 
tween those two regimes is naturally of the order of r s w 10 
when W is small. When W is larger, (AE) oc rf as indi- 
cated by the dashed lines of Fig. El Assuming that the 
correlated regime occurs at the ratios rf(W) below which 
(AE) does not depend on r s , and above which (AE) in- 
creases as rf, we have plotted in the inset of Fig. 1141 how 
rf(W) depends on W. The onset rf(W) of the correlation 
effects decays from the value « 10 where there is a level 
crossing without disorder towards a much smaller value 
when W — > 5. 

This onset rf(W) for W — 3... 5 can also be seen 
in the sample-to-sample distribution of the first energy 
excitation. When r s is weak, the first many body excita- 
tion corresponds to a single one body excitation above the 
Fermi energy. If the one body motion is diffusive, the one 
body spectrum is correlated, and the distribution P(s) of 
the energy spacing between consecutive levels is given by 
the Wigner surmise Pw(s). Therefore, the first many body 
spacing s — E\ — Eq/ (E\ — Eq) measured in units of its 
ensemble averaged value is given by the Wigner surmise: 



Pw(s) = — exp( — ) 



(21) 



as one can see in Fig.Elf° r r s < r f (W). When r s exceeds 
rf(W), the level repulsion becomes weaker. As shown in 
Fig- El P( s ) for l° w s is systematically larger than Pw(s) 
when r s = 5. To study this weakening of the spectral 
rigidity, we define a spectral parameter 



V(r s ) 



J a (P( S )-P w (s))ds 
f a (P P ( S )-P w ( S ))ds> 



(22) 
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Fig. 15. Distribution P(s) of the first energy excitation s = 
Ej - E / (Ej - Eo) for W = 5 and r s = 1.7 (•) and r s = 5 
(□). The solid line is the Wigner surmise Pw{s). 



value rf (W) than in the clean limit, yielding for the ex- 
citation energy AE an increase of its average value when 
r s increases and a change of its distribution. 

4.3 Breakdown of the Hartree-Fock approximation 

A last evidence proving that rf(W) is indeed the onset 
of the correlation effects is shown by Fig. El where the 
disorder average total GS persistent current (I x ) enclosing 
a flux <j) — 7r/2 is shown as a function of r s . Both the 
exact numerical value and the mean field value given by 
the Hartree-Fock approximation are given. As one can see, 
the change seen in the mean and in the distribution of the 
first energy excitation at rf(W) are correlated with the 
fact that the ground state cannot been described by the 
best possible Slater determinant (HF-approximation), due 
to correlation effects which are beyond a simple mean field 
approach. 




Fig. 16. Spectral parameter rj for the first energy spacing 
as a function of r s for W = 3 (?) and 5 (•), showing the 
sharp breakdown of the Wigner distribution at rf ~ 1.7. The 
inset shows that the level repulsion is totally suppressed when 
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Fig. 17. Ensemble average GS current (I x ) as a function of r s 
for W = 5. Exact values (filled symbols) and HF values (empty 
symbols) . 



where Pp(s) = exp— s is the Poisson distribution charac- 
terizing uncorrelated levels and a — 0.4729 is the value 
where Pw(a) — Pp(a). r\ = 1 when P(s) = Pp(s) and 
r] = when P(s) = Pw{s). Very remarkably, one can 
see in Fig. ^| that the first energy excitation is well de- 
scribed by the Wigner surmise up to a threshold consistent 
with rg (W) where the spectral rigidity suddenly becomes 
weaker. The curves rj(r s ) are given on a larger interval of 
values of r s in the inset of Fig. El an d one can see the two 
characteristic thresholds detected in earlier studies |14ll5j 
for W = 5: rf(W = 5)~2 where one has the breakdown 
of Wigner-Dyson rigidity and r]f (TF = 5) 10 where was 
located the onset of charge crystallization. 

In summary, when disorder yields one particle diffu- 
sion, the onset of correlation effects occurs at a weaker 



5 Conclusion 

Firstly, it is interesting to compare our results obtained 
on a 2d torus without edge to those obtained using an 
harmonic confinement. For instance, a Monte Carlo study 
|26| of a few electrons in an harmonic trap concludes that 
mesoscopic Wigner crystallization proceeds in two stages: 
(i) via radial ordering of electrons on shells and (ii) freez- 
ing of the intershcll rotation. This crystallization in two 
steps, with a particular intermediate behavior, could be 
attributed to the non uniform density characterizing the 
harmonic trap. One may argue that crystallization takes 
place in the low density edges before the large density 
bulk, such that this intermediate regime might be related 
to some interplay between edge and bulk effects. We find 
that mesoscopic Wigner crystallization takes also place in 
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two stages when the particles are confined on a 2d torus 
with a uniform density. 

Secondly, this raises the question to know if this in- 
termediate regime is a pure mesoscopic effect valid only 
in small systems, or the mesoscopic trace of the Andreev- 
Lifshitz supersolid, where unpaired fermions with reduced 
Fermi energy co-exist with a floppy solid. This might ex- 
plain a few recent experimental studies 7 8 9 of the metal- 
lic 2d hole gas in GaAs, where a not well identified metal- 
lic phase seems to coexist with a more usual Fermi liq- 
uid phase, responsible of usual weak localization behav- 
iors. This two-phase coexistence scenario may be simply 
explained by a supersolid phase, without having to nec- 
essarily take into account the disorder effects |57]. Upon 
completion of this manuscript, we received from Boris Spi- 
vak |28) a preprint where the existence of an intermediate 
phase between the Fermi liquid and the Wigner crystal is 
claimed to be a generic property of the 2d pure electron 
liquid in MOSFETs at zero temperature, and where the 
consequences for the experimental results obtained in 2d 
MOSFETs are also discussed. 

Thirdly, one cannot exclude that a supersolid regime 
is favored by the chosen geometry, because of the number 
of particles and underlying square lattice, but not favored 
at all in the continuous limit, which has not a square sym- 
metry but a spontaneously broken hexagonal symmetry. 
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